Tunnelling magnetoresistance anomalies of a Coulomb blockaded quantum dot 
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We consider quantum transport and tunneling magnetoresistance (TMR) through an interacting 
quantum dot in the Coulomb blockade regime, attached to ferromagnetic leads. We show that 
there exist two kinds of anomalies of TMR, which have different origin. One type, associated with 
the TMR sign change and appearing at conductance resonances, is of a single particle origin. The 
second type, inducing a pronounced increase of TMR value far beyond 100 %, is caused by electron 
correlations. It is manifested in-between Coulomb blockade conductance peaks. Both the types of 
anomalies are discussed for zero and finite bias and their robustness to the temperature increase is 
also demonstrated. The results are presented in the context of recent experiments on semiconductor 
quantum dots in which similar features of TMR have been observed. 

PACS numbers: 85.75.-d, 73.23.Hk, 73.63.-b 



I. INTRODUCTION 

Spin dependent tunnelling phenomena have attracted 
much scientific attention recently mostly due to the 
promising potential applications for magnetic sensors and 
magnetic random access memories. Early work of Mood- 
era et. al}, showing large, reproducible TMR effect of 
ferromagnetic tunnel junctions with AI2O3 spacer mea- 
sured at room temperature, was one of the first steps ini- 
tiating an avalanche of both theoretical and experimental 
investigations. The results were well understood within 
JuUiere's model^. Within this model the TMR effect is 
described in terms of the densities of states polarizations 
at Fermi energy of left (L) and right (R) ferromagnetic 
electrodes: TMR = 2PlPr/(1 - PlPr)- The electronic 
structure of the spacer is assumed to be featureless. It 
appeared that this model was too simple to understand 
the experimental results that further arose. The interface 
resonant states appearing as a result of energy structure 
mismatch of the ferromagnetic lead and insulating spacer 
have considerable influence of the TMR value^'**. In par- 
ticular, symmetry of these states can select spin-polarized 
bands in ferromagnetic electrode and enhance tunnelling 
through the insulating barrier. Moreover, a modification 
of the energy structure of the spacer and its bonding to 
the ferromagnetic leads can switch the polarization of the 
current and also change the sign of tunnelling magnetore- 
sistance. These features can be controlled by a proper 
modification of the spacer composition^. It was also 
shown that the defect states, present in the insulating 
spacer separating two ferromagnetic leads, can dramati- 
cally change the exchange interaction between the leads^ 
and affects TMR in the way not understood within Jul- 
liere's model. Recently it has also been derived that the 
scattering of electrons on nonmagnetic impurities present 
in the nanojunction barrier can cause the switching of 
TMR sigriS.. 

In recent years, due to a rapid development of nan- 
otechnology, a new kind of " spacer" became available to 
put between ferromagnetic leads. Importantly, the elec- 



tronic properties of these "spacers" are well controllable. 
These are semiconductor quantum dots^ (QDs). They 
have well defined discrete energy spectrum, whose posi- 
tion with respect to the chemical potential of the leads 
can easily be capacitatively tuned by electric field of a 
nearby gate. This in turn gives the possibility of elec- 
trical control of the TMR effect for such a device, which 
opens new possibilities of applications. Manipulation of 
the spin by electric field is one of the central issues of 
spintronics, which is regarded as a promising alternative 
for traditional charge-based electronics. Electron interac- 
tions inside the dots cause dramatic effects in their con- 
ductance, displaying quantum Coulomb blockade^ and 
Kondo effect^. Recently TMR measurements have been 
performed for such designed devices: InAs quantum dots 
coupled to nickel or cobalt electrodesiSiiiii^iH. These 
experiments show a rich TMR behavior, including the 
changes of the TMR sign and appearance of its maximum 
far exceeding 100 %. Spin transport and gate control of 
the tunnelling magnetoresistance has also been realized 
in carbon nanotube o^^'^^ . The variety of TMR anomalies 
observed in those systems is also ascribed to the discreet- 
ness of the nanotube energy structure. 

Tunnelling magnetoresistance oscillations caused by 
the classical Coulomb blockade at a small metallic droplet 
coupled to ferromagnetic electrodes had already been 
predicted long agc^^. The TMR value can also change 
its sign due to strong electron correlations inside the dot 
in Kondo regime, as was shown theoreticallyi^ii^. It is 
a result of the current enhancement by the Kondo reso- 
nance when the dot is coupled to the leads of antiparallel 
spin configuration. 

In the present paper we make an attempt to describe 
experimentally encountered anomalies of TMRiSi^iii^ for 
the InAs quantum dots in Coulomb blockade regime. We 
utilize the model of an interacting quantum dot with one 
level active in transport, coupled to ferromagnetic leads. 
We show that TMR sign switching is caused by the res- 
onances of the dot level with one of the leads chemical 
potential in presence of the large asymmetry of the dot- 
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leads coupling. We also predict that electron-electron 
interactions inside the dot have a decisive role in the 
formation the TMR maximum, which exceeds 100 % . 
Recently, within the similar model device, we also intro- 
duced a proposal of the electrical control of the spin po- 
larization of the current'''-. We will show that the correla- 
tion induced switching of the spin direction of the current 
is closely related to the TMR sign change at Coulomb 
blockade. 



II. THEORETICAL APPROACH 

The device is described by Anderson hamiltonian^*^ , 
where the dot takes the role of magnetic impurity and 
the (polarized) leads are analogues of host metal: 

H = edd^da + Un^Tia + ^ [^acl^^^d^ + h.c] 

+ ^ ] ^ka,(jC^a,<j^ka,(j (1) 

The first two terms describe the dot with the presence of 
Coulomb interactions U . The bare dot level is shifted 
by the gate voltage acting on the dot capacitatively: 
ed = — Vg, and its initial position for Vg = is as- 
sumed to coincide with Fermi level e^; = = 0. The 
third term describes the tunnelling between the dot and 
the leads, represented by the last term in Eq. The 
electron energy in the leads is spin-dependent, a =t,J., 
because the leads are assumed to be spin polarized. We 
neglect the spin dependence of the tunnelling matrix 
elements ta (a = R) which are rather dependent 
on the potential barrier between the dot and a given 
lead. Thus, the spin dependence of the QD level width 

(r<,/2) = (l/2)X;araa; = 27r|t„|2p„^ is caused by resonan^s'T^ =Tw and T^'^'' 



The retarded dot Green function G^(t — t') = ~i0{t — 
t'){da{t)dl{t') +dl{t')da{t)) is obtained by solving the 
set of equations of motion of the Green functions in the 
Hubbard I approximation^^. Within this approximation 
the two-particle Green functions describing spin-flip pro- 
cesses (generating Kondo effect) on the localized level are 
neglected. The Green functions that describe the normal 
scattering of band electrons on an impurity are approx- 
imated by decoupling of band electrons from impurity 
electrons. The Hubbard approximation is valid for large 
U/T ratio, when the Hubbard subbands are well sepa- 
rated in energy scale. For numerical calculations we as- 
sumed = 0.3C/, the other width components are cal- 
culated from Eqs. ^ and ^ for given lead polarization 
and asymmetry a. 

The Hubbard approximation is the simplest scheme 
which describes correlated electrons, placed on the 
approximation scale between Hartree-Fock approxima- 
tion for interacting but uncorrelated electrons, and the 
schemes for strongly correlated electrons, leading to 
Kondo physics. Thus, it is most suitable for the descrip- 
tion of a spin-degenerate QD level in the Coulomb block- 
ade regime of the lead-dot coupling, the limit realized in 
recent experiment a^°'^^'^^ . 

The Fourier-transformed expression for QD Green 
function with the spin a =1, i for given f3 — P or AP 
arrangement has the form: 



g;'^(^) = [ 

1 - (n^)^ 
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(4) 



the coupling to the leads with different spectral densi- 
ties Pa] 7^ Pal, which are assumed to be featureless and 
constant. 

Let us define the polarization of the quantity X, Px = 
[X] - X])/{X] + Xi). For the lead a it is: Pa = [pa] - 
Pai)/{pai + Pal), which cau be expressed by the spin- 
dependent QD widths: 



Eq. ([U has been written as the sum of two hubbard 

ed + U, whose spectral 



Pa — i^al - ^al)/ {^al + ^al)- 



(2) 



To calculate TMR we will consider parallel (P), Pr — Pl, 
and antiparallel (AP), Pr ~ —Pl, leads polarizations 
arrangement The asymmetry of the dot-leads coupling is 
described by a parameter. Because of validity of Eq. ([2]), 
the relations between dot level width components from 
left and right lead follow for (P) and (AP) configurations: 



i Ra ^OiL^g, 



(3) 



Tunnelling magnetoresistance is calculated from the for- 
mula: TMR = {0^ — Q^^)/Q^^ , where Q are appropri- 
ate conductances calculated for parallel and antiparallel 
configurations. 



weights are controlled by the dot level occupancy with the 
opposite spin a. This feature directly reflects Coulomb 
interactions between electrons with opposite spins. 

The spin components of the dot occupancy have been 
calculated selfconsistently from the set of coupled equa- 
tions: 

= j G<'^{cj,{n,f)du;, 
{n,f j Gf\uj,{n,)P)du:. (5) 

The "lesser" dot Green function C^'^ can be ex- 
pressed by the spectral density of the dot^^, p^{(-o) = 
-(l/7r)3G;./3(a;), G<'P{lo) = 2z^/>)pg(c.). Non- 
equilibrium distribution function / = [^LafL + 

^Rtrful/O^La- + ^Ra) ^ two-stcp profllc dcflucd by 
the chemical potential in the leads: /l/r. = T eV) 
and collapses into equilibrium Fermi-Dirac distribution 
function f = = in the limit of zero bias between 
the leads, eV 0. The current is calculated within Lan- 



3 




The conductance has exactly the same form as for non- 
interacting dot level of Green's function: = [oj — + 
i{r^/2)]~^ crossing Fermi level, — ep- Thus, these 
TMR features are of the single particle origin and can 
be described in the limit of non-interacting electrons. In 
Fig. Il]) also the TMR curve for non-interacting dot level, 
calculated for a = 0.1, is shown. The minimum at Vg = 
coincides with the one of interacting case for first hubbard 
level in resonance with Fermi level = ep. 

Taking into account the relations between the level 
widths, Eq. ([3]), it can be shown by straightforward calcu- 
lation that for = or e^/ = the spin components of 
conductance for parallel and antiparallcl configurations 
are (in units of /h): 



4a 



{1 + ay 



(10) 



FIG. 1: TMR dependence on gate voltage calculated for 
T = 0.01(7, Pl = 0.5 and zero bias for various asymmetry 
parameters :a = 0.1- solid, a = 1-dotted, a = 0.2 dash- 
dotted curve. The dashed line is for non-interacting dot and 
a = 0.1. 



dauer formalism from the relation^^: 



-^13 



L<T^ ^ Rcr 

In the limit of zero bias the conductance has the form: 
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III. 



dV 



-T. 



duj{- 



9/. 



-p/3 -p/3 



■pfH- (7) 



Rcr 



BEHAVIOR OF THE SYSTEM AT ZERO 
BIAS 



and 



±(1. 



(11) 



where also the relation between spin-dependent widths 
from Eq. Tli = -Fl|(Pl - 1)/(Pl + 1) has been 
applied. 

For symmetric coupling, a = 1, the conductance for 
(6) parallel configuration — Q"^ + reaches its maxi- 
mum value of 2e'^/h. It is reflected in the TMR curve 
in Fig. ([T]), which also displays a local maximum at the 
conductance resonances at Vg = and Vg — U . The 
conductances for antiparallcl configuration and a = 1 



are: Q 



AP 

T 



{Pl - I), which yields TMR 



P£/(l — Pl) > 0. Thus, for symmetric dot-leads cou- 
pling the TMR has positive sign at the conductance res- 
onances. The situation changes for asymmetric coupling, 
a < 1. In this case TMR at resonances takes the form: 



For all numerical results presented we have chosen the 
left lead polarization P^ = 0.5 and temperature T = 
O.OIC/. This range is typical for experimentsiSiiiii^; it 
gives the temperature of nimK for U — 15meV—. 

In Fig. ID) calculated TMR evolution with the change 
of gate voltage for different a asymmetry parameter is 
demonstrated. In order to understand various TMR 
anomalies shown in Fig. ([1]) it is instructive to analyze 
analytically TMR expression at T = 0. The conductance 
for a given spin and leads configuration f3 has the follow- 
ing form in this limit: 



7/3 







Ra 



r ed(ed + j7) 
Ud-(-C/(l-(n^)f) 



4^-'^ La 



(8) 



Ra) 



Consider the situation when or e^J hubbard level 
crosses Fermi level (at Vg — or Vg = U). In this case 
Eq. ^ 



takes the form: 



rP p/3 

^ La^ Ra 



h ifr'^ 

4 Lo 



(9) 



TMR : 



[^1(1 



(1 



^212 



(1 



HP|-l)[P|(l-a)2 + (l-fa)2 



(12) 

and for a -C 1 it changes the sign to negative: TMR — 
~2Pll{Pl -I- 1) < 0, as shown in Fig. 

Let us summarize the above discussed single particle 
TMR sign changes within a simple physical picture. For 
the perfect symmetric coupling, a — \, the transmis- 
sion through the dot in parallel configuration reaches 
the conductance quantum in the both spin channels, 
= Q^ ^ jh. This is caused by the perfect matching 
of the spectral densities of spin up and spin down of the 
left lead to the corresponding spectral densities of the 
right lead at the Fermi level. It also results that the spin 
dependent level widths due to the coupling to the left lead 
and to the right lead are equal, see Eq. ([3]). In antiparallcl 
configuration and a — 1 this is not the case, the numbers 
of states of spin up and spin down in the left lead at the 
Fermi level are different as compared to the right lead and 
also the corresponding level widths are not equal. Thus, 
the conductance in antiparallcl configuration is less than 
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conductance quantum, Q^^ = Q"^^ — (3/4)(e^//i) for 
Pl = 0.5, and TMR > 0. Note that for unpolarized 
leads, Pl — 0, the matching of the spectral densities is 
retained and the conductances Q:^^ — Q^^ = e^/h reach 
conductance quantum. When the asymmetry of the dot- 
leads coupling is increased, a < 1, the tunnelling between 
the dot and the right lead is reduced, which destroys 
the perfect matching of the widths in (P) configuration: 
rj^l 7^ r£| and r|^| 7^ r£|. it results in a gradual de- 
crease of the both spin conductance components in (P) 
configuration with the increase of a. More interesting 
situation takes place in (AP) configuration. The Q^^ 
conductance component, which describes the tunnelling 
of spin up excess electrons (Pl > 0) from the left lead 
via resonant dot state into minority up spin subband of 
the right lead {Pr < 0) decreases with the increase of the 
coupling asymmetry because the width which is less 
than r^^i" even for a = 1 is further decreased by lowering 
a. It increases the mismatch between the widths in this 
spin sector. Different relation is encountered between 
the widths in the (AP) spin down channel. Initially, for 
a = 1 we have T"^^ > Tf^^ and by the decrease of a 
the value of is lowered. When it reaches the value 
T"^^ — Tf^^ we have a perfect symmetric coupling in 
this channel. Thus, the decrease of the dot-leads cou- 
pling symmetry causes an increase of the coupling sym- 
metry in the spin down sector when the system in (AP) 
configuration. The component Gf^ increases and starts 
to dominate over other conductance components causing 
TMR < for small a. It can be checked from Eq. (fTTjl 
that Q'^^ reaches full transmission limit for a = 1/3 when 
Pl = 0.5. This enhancement of the one of the conduc- 
tance spin components in (AP) configuration can easily 
be generalized to other (AP) arrangements. For instance, 
if the polarization of the left lead were assumed to be 
negative, the Q^^ would be enhanced by the decrease of 
a. 

Similar as discussed TMR sign changes have been ob- 
served in quantum dot&iiii^. and in carbon nanotubes^^ 
coupled to ferromagnetic electrodes. They are also inter- 
preted in terms of the asymmetry of the coupling to the 
leads. 

Consider now the region in which the dot is in Coulomb 
blockade. In this case for ej, = —U/2 exactly one electron 
is present at the dot -I- {ni)^ = 1. There are con- 

siderable TMR anomalies near this point (see Fig. ([T])), 
TMR can change its sign or can be greatly enhanced 
exceeding 100 percent. These anomalies are caused by 
electron-electron interactions . The peculiar behavior of 
TMR in this region is caused by interplay of the two 
factors. The first is caused by the coupling (and its 
asymmetry) of the dot to the spin polarized leads. It 
determines the widths of the conductance peaks for par- 
allel and antiparallel configurations. The second factor 
are electron correlations which are the strongest in this 
region^^. They cause the spin components of the oc- 
cupancies to be close to one half at Coulomb blockade. 



The TMR anomalies here can be understood by analyz- 
ing the spin components of the conductance for given 
lead polarization arrangement, Eq. ([H]). For ^ —U/2 
and {rig)^ ~ 0.5 the denominator + U{1 — {ng)) is 
very small, causing to be small. Depending on the 
leads polarizations and the strength of the Coulomb in- 
teractions the components of the conductances defining 
TMR get their minimal values for different positions of 
the dot level. The same mechanism causes sudden cur- 
rent polarization switching in the region of —U /2 as 
discussed in^^. Indeed, it is shown in Fig. ^ that for (P) 
and (AP) configuration the value of conductance polar- 
ization Pg performs a rapid oscillation and additionally 
it changes sign for (AP) configuration . 

TMR value in this region is very sensitive to the change 
of the conductance in (AP) configuration. Consider first 
the perfect symmetric coupling, a ~ 1. For P^ = —Pl 
the total widths of the dot level are equal F^ — F| be- 
cause of relation Tl^ — and Tli = F^i, Eq. 
Thus, in the symmetric (AP) arrangement the dot be- 
haves as if it were coupled to unpolarized leads. In such 
a case spin components of the conductance are equal, 
— Qy^^, in the whole range of gate voltages and also 
the occupancies (n|)'^^ = (n^)^^ are equal. The QD 
occupancy curve displays a plateau at (n^) = (n|) ~ 0.5 
due to Coulomb blockadei^. Moreover, for particle- 
hole symmetric case, = the occupancies are 
{n^)AP = {n{)^P = 0.5, giving gf^ = Qf^ = 0. This 
feature causes an infinite TMR value for a = 1 at T = 0. 
For asymmetric coupling, a < 1 and finite temperature, 
the relation Q:^^ — Q^^ is still fulfilled at the gate volt- 
age of TMR maximum, but the conductances have small 
finite value. From the condition of Q^^ = Q-^^ follows 

also the equality (n-^)"^^ — (ni)^^ at gate voltage of 
TMR maximum, independently on a. It can easily be 
derived utilizing Eq. ([5]). At Coulomb blockade the first 
term in the denominator is much larger than the second 
one and the relations Eq. ^ for (AP) configuration have 
to be used. Thus, at TMR maximum the conductance in 
(AP) configuration is unpolarized (the conductance po- 
larization Pg^ = as shown in Fig. ([2])) and also the 
dot occupancy polarization is zero. 

Strong enhancement of TMR at classical Coulomb 
blockade has also been predicted for ferromagnetic dou- 
ble tunnel junctions^ as a result of cotunneling of elec- 
trons through metallic island. It has also been observed 
experimentally^-?;^?,. Our result provides an explanation 
of corresponding TMR maximum at quantum Coulomb 
blockade for semiconductor quantum dot. It is in rela- 
tion with the recent TMR measurements for InAs quan- 
tum dot coupled to the Ni leadsi^, where the TMR en- 
hancement above 300 % has been observed at Coulomb 
blockade. 

At the gate voltages for which TMR — following 
equalities apply Qf^^ = G^ and — . It im- 

plies that the conductance spin polarizations are opposite 
for parallel and antiparallel configuration: Pg ~ —Pg^ 
(compare corresponding curves in Fig. 
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Let us now discuss the TMR sign change encountered 
at Coulomb blockade. It is shown in Fig. ^ that this 
TMR minimum coincides with the sharp minima of con- 
ductance polarizations, moreover Pg^ < here. The 
TMR minimum is caused by sharp polarization switch- 
ing of the conductance in (AP) configuration. It is in 
contrast to (P) configuration, for which the conductance 
polarization does not change its sign. In (P) configura- 
tion the dot is coupled to the leads both having excess 
of electrons of the same spin (in our case spin up) and 
conductance of this spin dominates in the whole range of 
gate voltages, Pg > (Fig. In (AP) configuration, 
on the average, there is no such an excess of electrons of 
particular spin coming from the leads. Thus, Q^^ and 
Gf^ are comparable in magnitude and the interactions 
between spin up and spin down electrons cause more dra- 
matic changes in the conductance polarization. Apart 
from electron-electron interactions, which cause the con- 
ductance polarization switching in (AP) configuration, 
the mechanism of enhancement of the conductance com- 
ponent Gf^ for small a works also here, similarly as dis- 
cussed for single particle origin TMR sign change. The 
transmission in the spin down channel of (AP) configura- 
tion dominates here over other conductance components 
and causes TMR < 0. There is no unique condition for 
negative TMR minimum to appear at Coulomb blockade 
and it depends on the value of initial leads polarization. 
For higher Pa, the conductance in (P) configuration will 
dominate {Pg curve then is shifted upward) and TMR 
will not change sign in spite of conductance polarization 
switching in (AP) configuration (not shown). The polar- 
ization P^^ behaves similarly as in the case of configu- 
ration with one lead polarized^^. The TMR sign change 
here also fades out rapidly with increase of temperature, 
Fig. ([31), because the Pg^ is very sensitive to the temper- 
ature change in this region. At higher temperatures Pg^ 
switching decays quickly. This sensitivity in the region 
of Vg ~ U/2 has been shown in^^ for one of the leads 
polarized. 

It could be counterintuitive that the conductance po- 
larization Pg reaches values larger than the leads polar- 
izations Pa = 0.5 for \Vg\ <^C ep, as displayed in Fig. 
We demonstrate that it is the case. In this limit the dot 
level is placed far from Fermi level, thus the dot is un- 
occupied or fully occupied. Consider the unoccupied dot 
for configuration which is realized for large negative 
gate voltage. Inserting {rig)^ — into Eq. ([5]) and noting 
that (erf)^ 3> (rcr/2)^ the conductance polarization can 
be written as: 



/3 



p/3 p 



r/3 r/3 



p/3 p/3 



p/3 p/9 
Ll'^ Rl 



(13) 



For (5 = P the relations of Eq. ([3]) are further utilized to 
give: 



r2 



r2 



2Pr 



Pi 



(14) 




FIG. 2: TMR (solid curve) vs. gate voltage calculated for 
a = 0.1 at r = 0.0117 and zero bias. The corresponding con- 
ductance polarizations Pg (dashed curve) and Pg^ (dotted 

r>A" 

Q 




FIG. 3: Temperature dependence of the zero bias TMR cal- 
culated for Q = 0.1 and Pl = 0.5: T = 0-dotted, T = 0.0R7- 
solid, T = 0.02t/-dashed and T = 0.03;7-dash-dotted curve. 



The last r.h.s expression has been obtained by substitut- 
ing Vli calculated from Eq. Note that the result is 
independent on a. Thus, for Pl — Pr — 0.5 the con- 
ductance polarization Pg = 0.8. The same expression, 
Eq. ([M]) . is obtained for fully occupied dot (rig)^ — 1 and 
noting that (e^ + f/)^ :$> For antiparallel con- 

figuration, the conductance polarization Pg^ = in the 
limit of iVgl <C Ci?, as shown in Fig. It can straight- 
forwardly be derived from Eq. (jl3p for f3 = AP when 
Eq. ([3]) is utilized: Note that TMR in the discussed limit 
of I Ed I <C ep reaches value 2/3 as predicted by JuUiere's 
model^ for two polarized leads separated by featureless 
tunnel barrier. 

Temperature dependence of the TMR for asymmetric 
coupling a — 0.1 is shown in Fig. (3). TMR anoma- 
lies in the regions of Vg and Vg ^ U caused by 
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FIG. 4: TMR dependence on gate voltage (solid curve) for 
finite bias eV = O.IU calculated at T = O.OlU, Pl = 0.5 and 
asymmetric coupling to the leads a = 0.1. The correspond- 
ing differential conductances (dJ/dV)^ (dotted curve) and 
(dJ/dV)^'' (dashed curve) are also shown. 



the resonances of QD hubbard levels and e^/ with 
the Fermi level are robust to the increase of tempera- 
ture. The anomalies due to electron correlations situ- 
ated in the range of Coulomb blockade Vg ~ J7/2 in 
turn, are sensitive to these changes. It is caused by 
a temperature increase of the conductance spin compo- 
nents at the Coulomb blockade valley. However, the pro- 
nounced TMR maximum reaching 100 percent survives 
for T = OmU (~ 0.5A' for U = IbmeV^). Contrary, 
the negative TMR minimum rapidly disappears at higher 
temperature, as pointed out previously. 



IV. BEHAVIOR OF THE SYSTEM AT FINITE 
BIAS 

If the finite bias is applied to the system, an electron 
transport through excited states of the QD can be acti- 
vated. In the Hubbard approximation, used for descrip- 
tion of the device, these processes are not taken into ac- 
count. 

We consider two limits of the value of the bias applied 
as compared to the width of the QD level: small bias 
limit: \eV\ < and large bias limit when \eV\ > F^. 
In the large bias limit the chemical potentials of the 
leads are well separated in energy scale. It implies also 
a good separation of the differential conductance reso- 
nance peaks which appear when any of the hubbard QD 
levels crosses given chemical potential. Thus, the max- 
ima and minima of TMR are also well separated. In the 
small bias limit the conductance resonances overlap each 
other which causes splitting of TMR minima and dimin- 
ishing of TMR maximum at Coulomb blockade. In both 
the limits, a correspondence can be found between TMR 
features with those appearing at zero bias. 



A. Small bias regime 

The TMR behavior for eV = Q.1U is shown in Fig. © 
for asymmetric coupling to the leads a = 0.1 and tem- 
perature T — O.OIC/. One notices that the pronounced 
TMR minima present for = and Vg = U are split 
when a finite bias is applied. For positive bias the left 
lead chemical potential /iL is shifted upwards by eV and 
right lead chemical potential is shifted downwards by 
eV on energy scale. When gate voltage increases from 
negative values, it shifts the dot level from empty state 
regime towards Fermi bathes inside the leads. At first 
the level comes into resonance with and then with 
fj,R chemical potential. It causes the appearance of two 
peaks in conductance with the distance of doubled bias 
value between them (shown in Fig. (|4|)) and consequently 
two minima of TMR are produced, labelled by (1) and 
(2). Similarly, in the range ofVg ^ U the second hubbard 
level e^J comes into resonance with fi^ and then with iiji 
and minima (3) and (4) appear. 

The maximum of TMR at Vg ~ U/2, distinct in 
case of zero bias, is diminished when the bias is finite. 
This maximum has been associated with the equality of 
Q:^^ = Q^^ and subsequently (n|)^^ = (?^x)^^ dis- 
cussed for eV = 0. For finite bias this relation is still ful- 
filled. It can be checked by analyzing spin components of 
the differential conductance in Coulomb blockade region. 
However, the conductance for antiparallel configuration 
is enhanced by finite bias and the splitting and shifting of 
the conductance peaks also appears. It causes a gradual 
diminishing of TMR maximum. The negative TMR min- 
imum present for zero bias and associated with the rapid 
change of the conductance polarization Pg^, has disap- 
peared at finite bias. It has been shown in^^ that this 
sudden conductance polarization switching also rapidly 
diminishes when finite bias is applied. 

B. Large bias regime 

Consider the case of eV = U for which the left (right) 
chemical potential fj,L (M/j)is shifted upwards (down- 
wards) by U on energy scale. The TMR vs. gate volt- 
age for such a bias and a = 0.1 is presented in Fig. ([5]) 
along with corresponding differential conductances. Let 
us discuss various anomalies appearing in the TMR curve 
labelled by arabic numbers. The TMR minima placed at 
the differential conductance resonances have single par- 
ticle origin. For instance, in the point (1) at Vg = —U 
the first hubbard level coincides with /i^ and the sec- 
ond e^/ lying above is empty. At this point TMR has 
a (negative) minimum. This kind of minimum has ap- 
peared already for zero and small bias each time when 
the QD hubbard level crossed chemical potential of the 
leads: for eV = it corresponds to the minimum at 
Vg = (see Fig. ([Ij), and to the minimum (1) for small 
bias (Fig. (|3])). Similar correspondence is for minimum 
(5) at = 2U where the second hubbard level e^/ co- 
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FIG. 5: TMR dependence on gate voltage (solid curve) for 
finite bias eV = U calculated at T = O.Off/, Pl = 0.5 and 
asymmetric coupling to the leads a — 0.1. The corresponding 
differential conductances (dJ/dV)^ (dash-dotted curve) and 
{dJ/dV)^'' (dashed curve) are also shown. The dotted TMR 
curve is for symmetric coupling a = 1. 



1.0 




-0.4 I 

-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

eV/U 



FIG. 6: TMR dependence on applied bias for the set gate 
voltage Vg = Q and Pl = 0.5, calculated at T = O.Olf/, and 
a = O.I-solid curve and a = 1-dotted curve; dashed curve- 
TMR for a = 1 and T = 0.06(7. 



incides with and is fully occupied. This minimum 
is analogous to the minimum at Vg = U for zero bias 
and minimum (4) for small bias in Fig. ([4]). The min- 
imum (4) at = where e^/ in resonance with /i^ 
corresponds to the minimum (3) for small bias. One 
can also note series of maxima and minima appearing 
in-between conductance resonances; for instance mini- 
mum (2) and maximum (3). These anomalies are caused 
by electron correlations and correspond to similar fea- 
tures of TMR at zero bias (Fig. H])) at Vg ~ U/2. For 
large bias, eV 3> F^, when the conductance resonances 
are well separated in energy scale, the Coulomb block- 
ade TMR anomalies can be identified and their corre- 
spondence with zero bias TMR features can also be es- 
tablished. In the small bias regime. Fig. they are 
considerably diminished. 

The TMR minima of single particle origin, for instance 
those numbered by (1), (4) and (5), transform into local 
maxima for symmetric dot-leads coupling, a = 1 (dot- 
ted curve). Similar feature was present for zero bias 
case. Fig. Thus, by experimentally tuning the cou- 
pling asymmetry one can resolve the mechanism causing 
a given TMR anomaly. Note that the TMR curve for 
a = 1 is symmetric with respect to Vg — U/2. At this 
point the hubbard levels lay in the middle of the trans- 
port window of the width 2U. For instance, when Vg is 
set to zero (to U), the e^J (e^) comes into resonance with 
Mi (ma) ^^^d the other hubbard level is shifted to the 
center of the transport window. It gives the same TMR 
feature for ±Vg. 

Consider now the bias dependence of TMR for the 
set gate voltage. In Fig. ([6]) the TMR bias dependen- 
cies for a = 0.1 and a = 1 calculated for = are 
shown. Again, the TMR minima labelled by (1), (2) 
and (3), present for large coupling asymmetry a = 0.1, 



can be identified as of single particle origin. For min- 
imum (2) at zero bias is in resonance with both 
and /i/j. For minimum (1) at eV = —U the e^^ hub- 
bard level is in resonance with fiL chemical potential 
and minimum (1) at eV ~ U the e^^ hubbard level is 
in resonance with /ifl chemical potential. These three 
minima disappear for symmetric coupling a — 1 and the 
local maxima develop instead as discussed before. The 
minima and maxima in the regions of eV ~ ±C/ /2 cor- 
respond to electron correlations induced TMR anomalies 
discussed previously for zero bias, Fig. |T]). These anoma- 
lies appear close to particle-hole symmetric case, when 
the Fermi level is situated in-between dot hubbard lev- 
els. In the present case, at eV = —U/2 the chemical 
potential is placed in-between hubbard levels in the se- 
quence: fXR < f-d < ^'L < f^d^ ■ At ey = U/2 the is in- 
between hubbard levels in sequence: /i^ < < < e^^ 
in energy scale. From comparison of the TMR curves 
for a = 1 and T = Q.OlU and at a higher temperature 
T — 0.06C/ one notes that the minima caused by electron 
correlations are much more sensitive to the temperature 
increase than maxima, similarly as in the case of zero 
bias. Fig. ([3]). Thus, by the increase of temperature one 
can distinguish the TMR minima of single particle origin 
which are robust to the temperature change from those 
caused by electron interactions, sensitive to temperature. 

For symmetric coupling to the leads the TMR curve 
becomes symmetric with respect to eV ~ 0. It is un- 
derstood when one notices that the change of the bias 
direction is equivalent to the simultaneous exchange of 
Ml ^ Mi?, ^-iid e^^. It gives for symmetric coupling, 
a = 1, the same value of TMR for zteF. 
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FIG. 7: Dot magnetic moment vs. gate voltage for parallel 
(a) and antiparallel (b) configuration for various asymmetry 
a. Panel (a) and (b): q = 0.1 solid line, a = 0.5-dashed. 
Dotted line in panel (a) for a = 1 and in panel (b) for a — 0.9 
(jnQD = for Of = 1 in antiparallel configuration). Calculated 
for zero bias at T = O.Olt/, Pl = 0.5. 



V. INFLUENCE OF EFFECTIVE MAGNETIC 
FIELDS 

The dot attached to spin-polarized leads is influenced 
by two effective magnetic fields, which have different ori- 
gin, i) The hopping of electrons from the leads, which 
have an excess of one of the spin components, produces 
an occupancy polarization of the dot. Thus, the dot ac- 
quires magnetic moment TOq^ = {{n-^)^ — (n|)^)|/is|, 
where /is — g^sSz is the magnetic moment associated 
with spin-dependent (but degenerate) sub-levels of td 
and Sz = ±1/2. This effect can be regarded as if the 
dot were under influence of an external magnetic field 
and the Zeeman splitting of the dot energy level is pro- 
duced. However, the value of the resultant dot's mag- 
netic moment at given temperature is regulated by the 
relative position of the dot level with respect to the leads 
chemical potentials and electron-electron interactions. It 
has non-monotonic behavior as shown in Fig. ([7]). It 
is in contrast to the usual Zeeman splitting of the level 
by external field H^xt, for which the magnetic moment 
m = tanh(\iis\Hext/kBT) increases monotonically with 
increase of the field. Let us discuss the general features 
of the dot magnetic moment for parallel and anti-parallel 



configuration. Fig. Because Pl > the spin sub- 
levels of each dot hubbard level e^^ and e^^ (7 = /, //) 
have different widths Fj > Fj^. For gate voltage large and 
negative the e^J is empty and both the and e^^ are 
barely populated; the magnetic moment is small. The 
value of mQD is positive because F| > F| and is be- 
ing populated earlier than e^j^. When Vg increases, mqu 
also increases reaching the maximum. Further increase 
of Vg causes faster filling of sub-level because it is 
sharper than e^^. Thus, niQD decreases towards zero at 
= 0, where is in resonance with Fermi energy ep- 
Further shift of causes further increase of (ni) until 
a minimum of mgo is reached. Note, that the sequence 
of the filling of the first hubbard level e^, caused by in- 
equality of the Fj and F|, is additionally enhanced by 
electron-electron interactions, manifested by the spectral 
weight dependence '--^ (1 — (rig)) of e^, Eq. After 
reaching of toqz? minimum, the second hubbard level e^^ 
starts to be filled and toqz? increases towards zero due 
to the same mechanism (F^ > F|) as in the region of 
Vg ^ 0. The niQD = at Coulomb blockade corresponds 
to the particle- hole symmetry case of = —U/2 and 
(n|) = (ni) = 0.5 for unpolarized leads. Further increase 
of Vg causes a maximum of mg/) to appear, followed by 
fTT-QD = 0, when ej/ is in resonance with ep for Vg = U, 
and then minimum. Note that the mgo maximum is 
weakened, as compared to one at Vg ~ 0, and the mini- 
mum is enhanced. It is due to electron interactions: the 
spectral weight of the second hubbard level e^J is ^ (ng-) 
(see Eq. (U))). It also causes that the dot magnetic mo- 
ment remains negative for large positive Vg (compare to 
wqd > for large negative Vg). 

An increase of the a symmetry of the dot-leads cou- 
pling causes a gradual decrease of the dot magnetic mo- 
ment for (AP) configuration, panel (b). For symmetric 
case of a = 1 the magnetic moment is zero in whole 
range of gate voltages as has been pointed out in the dis- 
cussion of TMR maximum at zero bias. For a < 1 the 
mQo = line collapses into one Vg point at which the 
TMR maximum appears. 

ii) The second field is an effective magnetic field, Heff, 
produced by magnetized electrodes. Its value is not de- 
pendent on gate voltage but rather on the relative leads 
polarizations. While TMR measurement is performed, 
the external magnetic field is applied to the system, 
which changes relative leads polarizations. Due to the 
intentionally different shape anisotropy of the left and 
right lead, the magnetization of each lead responds dif- 
ferently to the external field. It enables anti-parallel lead 
polarization for small field and parallel configuration for 
higher field. The dot itself is in turn subjected to an effec- 
tive magnetic field produced by the magnetic electrodes 
nearby. This field is the largest for parallel configuration. 

For InAs quantum dots gyromagnetic ratio is negative, 
g < 0^^, thus the spin moments of e^i and e^x point in the 
same direction as the corresponding magnetic moments. 
The dot level subjected to the field Heff produced by 
the polarized leads (of direction opposite to the leads 
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spin polarization) is Zeeman split:ed| = — A and e^t = 
Erf + A, A = Ig^sHef f The Zeeman energy splitting 
for I5I = 3.8i^ gives 2A = 0.22meV/T {= O.OlSt/ per 
tesla) . We considered the influence of H^f / on the TMR 
features, assuming H^ff of the order of ±0.ir, taken 

from expcrimentsiSiiii. For antiparallel configuration the 
cancellation of the fields originating from the leads has 
been assumed H^^^ = 0. We have found that the possible 
Zeeman splitting of the dot level has negligible effect on 
TMR in this field range. 

VI. CONCLUSIONS 

We have discussed the TMR anomalies encountered for 
the quantum dot coupled to spin-polarized leads in the 
regime of Coulomb blockade. We have shown that there 
are two kinds of such anomalies. One kind has the single 
particle origin and can be interpreted in the frame of non 
interacting electrons model. The second kind of anoma- 
lies is caused by electron interactions. The TMR minima 
(and its sign change) of single particle origin appear at 
the conductance resonances for asymmetric dot-lead cou- 
pling. They are robust to the temperature increase and 
gradually transform into local maxima when the sym- 
metry of the dot-lead coupling increases. The anomalies 
associated with electron interactions appear at Coulomb 
blockade, in-between conductance resonances. The TMR 



maximum at Coulomb blockade, far exceeding 100 % is of 
this origin. It appears when in antiparallel configuration 
both the current the dot occupancy are spin unpolar- 
ized. This maximum survives at typical temperatures of 
experiment. We also predict the TMR sign change at 
Coulomb blockade. It appears due to the rapid polariza- 
tion switching of the current in (AP) configuration and 
the enhancement of the conductance in one of the (AP) 
spin channels by the dot-leads coupling asymmetry. It 
is very sensitive to the increase of temperature and de- 
pends on the initial polarization of the current coming 
from the leads. We have shown that the nature of the 
discussed anomalies can be experimentally resolved by 
the change of the dot-leads coupling asymmetry and/or 
temperature. Finally, we have analyzed the dot polariza- 
tion, as induced by the coupling to the polarized leads, 
and shown that it also depends on electron interactions 
present inside the dot. The estimated Zeeman field split- 
ting, produced by the leads, has negligible effect on TMR 
for experimental range of the fields. 
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